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ABSTRACT 

We present an XMM-Newton observation of the radio galaxy 4C 23.56 at z = 2.48 
which reveals extended X-ray emission coincident with the radio lobes spanning ~ 0.5 
Mpc. These are the largest X-ray-bright lobes known at z > 2. Under the assump- 
tion that these features result from inverse- Compton scattering of cosmic microwave 
background photons by relativistic electrons in the radio source lobes, the measured 
luminosity of Lo.s-gkeV = 7.5 x 10^^ erg s~^ implies a minimum energy stored within 
the lobes of ~ 10^^ erg. We discuss the potential of the large-scale radio/X-ray lobes to 
influence evolution of the host galaxy and proto-cluster environment at high redshift. 

Key words: galaxies: individual: 4C 23.56 — galaxies: high-redshift — galaxies: 
evolution — X-rays: galaxies: clusters — radio continuum; galaxies 



1 Introduction 

The galaxies hosting powerful radio-loud active galactic nu- 
clei (AGN) at high redshift are believed to evolve into the 
most massive galaxies in the present-day Universe, which 
reside locally in rich cluster environments (e.g.. Best et al., 
1998). High-redshift radio galaxies (HzRGs) are therefore 
popular targets in studies of the co-evolution of massive 
galaxies, their central AGN, and their environment on larger 
scales. 

In recent years, observations with Chandra and XMM- 
Newton have detected extended X-ray emission in an in- 
creasing sample of HzRGs (Fabian et al. 2003, 2003b; Carilli 
et al. 2002; Yuan et al. 2003; Scharf et al. 2003; Belsole et 
al. 2004; Overzier et al. 2005; Erlund et al. 2006; Blundell et 
al. 2006). Generally associated with the radio structure, the 
extended X-ray emission is believed to arise in most cases 
partly or wholly from inverse- Compton (IC) scattering of 
cosmic microwave background (CMB) photons by relativis- 
tic electrons in the jets and/or lobes of the radio galaxy. As 
noted by Schwartz (2002), the steep (1 + dependence 
of the density of CMB photons on redshift suggests this 
process will be increasingly significant toward higher red- 
shifts. Based on extrapolation from the radio source number 
counts, Celotti & Fabian (2004) estimate that at redshift s 
above z>l the majority of luminous {Lx > 10^^ erg s~^) 
extended X-ray sources may result from IC scattering of the 
CMB in HzRGs. 

This new class of X-ray source may provide crucial in- 
sight into the eflFect of AGN activity on the evolution of mas- 
sive galaxies and their larger scale environment. The iden- 
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tification of bubble-like features in low-and intermediate- 
redshift clusters indicates AGN jets can inject a significant 
amount of energy into their environs (e.g., McNamara et 
al. 2000; Fabian et al. 2006). Similar processes have been 
suspected to exert significant influence on the evolution of 
structure at high redshift, e.g., through the regulation of gas 
accretion and star formation in massive galaxies (e.g., Silk 
& Rees, 1998) and as a source of non- gravitational heating 
within forming cluster cores (e.g., Bower et al., 2001). How- 
ever, the mechanisms and magnitude of AGN feedback at 
high redshift are unknown. 

The detection of IC-scattered X-ray emission in HzRGs 
has begun to clarify the evolution and energetics of high- 
redshift radio lobes. As energy losses from the more ener- 
getic electron population responsible for the radio emission 
are more rapid than those of the less energetic population 
involved in IC scattering, the X-ray lobe emission can out- 
last the radio emission and provide a measure of past radio 
activity (e.g., Fabian et al. 2003; Blundefl et al. 2006). X- 
ray measurements have allowed estimates of the minimum 
Lorentz factors, 7min, in radio hotspots (Blundell et al., 
2006) and the total energy contained within X-ray lumi- 
nous radio lobes (Erlund et al. 2006) at z ^ 2. Scharf et al. 
(2003), noting the relative positions of X-ray and Lya flux 
in 4C 41.17 at z = 3.8, suggest that the IC-scattered X-ray 
photons could themselves give rise to a secondary feedback 
mechanism by photoionising the surrounding medium. Both 
the direct mechanical feedback from the jets and this second 
radiative mechanism could influence the structure and en- 
ergetics of forming structures on galactic and cluster scales 
and, interestingly, would preferentially aflFect the most mas- 
sive structures at high redshift. 

In this Letter we report the detection of extended X- 
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Figure 1. a: 0.5-8 keV EPIC image of 4C 23.56. b: Smoothed 2-8 keV EPIC image of 4C 23.56, with the region used for spectral 
analysis of the X-ray point source overlaid. The darker line seen through the NW corner is a chip gap in the pn camera, c: Smoothed 
0.5-2 keV EPIC image of 4C 23.56, with regions used for spectral analysis of the extended emission overlaid, d: Difference image obtained 
by subtracting a smoothed 0.5-2 keV band image from a smoothed 2-8 keV band image. Soft flux appears light (yellow) and hard flux 
appears dark (blue). Labels indicate the Knopp & Chambers (1997) 6 cm radio positions of the 4C 23.56 core (C) and NE (B) and SW 
(A) lobes. Images are 2' on a side and orientated so that North is up and East is to the left. 
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ray emission associated with the radio galaxy 40 23.56 {z = 
2.483; Rottgering et al., 1997) in XMM-Newton observations 
obtained as part of a larger study of the X-ray properties of 
HzRGs. We begin with a summary of the mult i wavelength 
properties of 40 23.56 in §2, present the X-ray observations 
and analysis in §3 and briefly discuss our findings in §4. 
Throughout this letter we assume a AODM cosmology with 
^m=0.3, Qa — 0.7, and Ho — 70 km s~^ kpc~^. Unless 
specified otherwise. X-ray fluxes are given in the observed- 
frame 0.5-8 keV band. X-ray luminosities are corrected for 
Galactic absorption and given in the rest-frame 0.5-8 keV 
band, and errors are quoted at the 90% confidence level. 

2 4C 23.56 

40 23.56 is an ultra-steep spectrum, Fanar off- Riley class II 
radio galaxy which is among the largest known at high red- 
shift. The core and radio lobes lie along an axis with a po- 
sition angle (P.A.) of 52° E of N and an angular size of 
roughly 0.5 Mpc (Chambers et al., 1996a). Multiwavelength 
observations indicate a heavily obscured AGN in a dusty 
host galaxy (Ohambers et al., 1996b; Knopp & Chambers, 
1997, Archibald et al., 2001). An enormous Ly-a emission- 
Une nebula comprised of two opposing ionization cones with 
opening angles of ^ 90° covers an area of ^ 5000 kpc^ paral- 
lel to the radio axis. The near-symmetry of the emission-line 
morphology indicates the radio axis lies nearly in the plane 
of the sky, while the relative brightness of the southwestern 
lobe implies it is somewhat nearer to us. Moderate excess 
counts of both Lya emitters and faint, red galaxies in the 
vicinity of 4C 23.56 suggest it may reside in an over-dense 
cluster or proto-cluster environment (Knopp & Chambers, 
1997). 

3 X-ray observations 

We obtained 32 ks XMM-Newton observations of the 4C 
23.56 field centred at 21^07^12.70", 23°3l'23.o" on 29 
November 2004 (ObsID 0204400201). The data were taken in 
Full Frame mode through the Thin filter and reduced using 
version 6.0.0 of the XMM Science Analysis System (SAS). 
The observations were minimally affected by periods of back- 
ground flaring, which we identified in 7 - 15 keV lightcurves 
and removed. The remaining useful exposures were 27.7 ks 
in the pn camera and 30.8 ks in the MOS cameras. We ap- 
plied the standard recommended PATTERN, FLAG, and 



energy filters to the data to produce event lists suitable for 
imaging and spectral analysis. 

3.1 Imaging analysis 

We produced images for each of the EPIC cameras in the 
full (0.5-8 keV), soft (0.5-2 keV), and hard (2-8 keV) bands 
by binning the filtered events into 2" pixels. We constructed 
a combined EPIC image by summing the individual instru- 
ment images. To improve visibility, we also created smoothed 
versions of all images by convolving the counts images with a 
Gaussian kernel roughly the size of the XMM PSF (a = 6'')- 

Figure la shows a 2^ x 2^ region of the unsmoothed, 
full-band, EPIC image centred on the position of 4C 23.56. 
The radio galaxy is clearly detected in the X-ray and ap- 
pears slightly extended toward the NE and SW. The ex- 
tended emission is more easily identified in the smoothed 
images of the same region. While the hard-band image (Fig- 
ure lb) shows only a point-source at the position of the radio 
galaxy, the soft-band emission (Figure Ic) is unambiguously 
extended along the NE-SW axis. An unrelated soft X-ray 
point source, visible in the unsmoothed counts image, is also 
seen at the southern edge of the Figure. 

To quantify the significance of the extended emission 
detection, we extracted hard and soft band counts within 
the pie shaped regions centred on the X-ray point source 
shown in Figure Ic. Two of the sectors are aligned with the 
axis of the extended emission while the other four serve as 
background regions. To the extent that the XMM PSF is 
spherically symmetric, these regions should contain roughly 
the same amount of contamination from the central point 
source. Across all three detectors, we find 144 soft-band 
counts within the two source regions and 156 counts within 
the four background regions. Correcting for area, we expect 
86 background counts within the source region, implying 
an excess detection within the source regions of 58 counts, 
which is significant at the 8.6cr level. The net source counts 
are split roughly evenly between the NE and SW regions. 
Considering the pn and the MOS counts separately, we find 
significant soft-band detections of the extended flux at the 
7.5cr and 4.8cr levels, respectively. We find no statistically 
significant excess of hard-band counts. 

To maximise the contrast of the faint, soft, extended 
flux, we subtracted the smoothed soft-band image from the 
smoothed hard-band image to produce the difference image 
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Figure 2. Top: X-ray spectra of the core of 4C 23.56 from the 
pn (upper points, black) and MOS (lower points, green and red) 
cameras, with the best-fitting model (histograms) described in 
the text overlaid. Bottom: Fit residuals. 



Figure 3. Top: X-ray spectra of the 4C 23.56 extended emis- 
sion from the pn (upper points, black) and MOS (lower points, 
green and red) cameras, with the best-fitting model (histograms) 
described in the text overlaid. Bottom: Fit residuals. 









— I-I — 1 — 1 — 1 — 1—1 










- Il i 1 1 1 ' 1 ' ^ 



channel energy (keV) 




channel energy (keV) 



shown in Figure Id. As both the central point source and the 
background have harder X-ray spectra than the extended 
emission, this method of displaying the data effectively de- 
emphasizes these sources of contamination, allowing the two 
lobes of soft emission to stand out distinctly on either side of 
the X-ray core. We stress that the morphology of this differ- 
ence image is not directly indicative of the surface brightness 
of the extended emission and must be interpreted with cau- 
tion. For example, the deficit of soft emission around the 
point source does not indicate the extended component is 
fainter near to the radio galaxy, rather that the contribution 
of hard flux in the wings of the point source is dominant over 
the soft emission at these radii. 

We indicate in Figure Id the 6 cm radio positions of 
the core of 4C 23.56 and the peaks of its NE and SW radio 
lobes (Knopp & Chambers, 1997). The peak of the X-ray 
point source is coincident with the radio position of the core 
and the bulk of the extended X-ray emission is aligned with 
the radio lobes at P.A.^ 52°. The peak of the excess soft 
flux coincides with the NE radio lobe at a distance of ^ 20" 
from the radio core 160 kpc z — 2.483). In the SW, the 
excess soft flux appears to peak interior to the radio lobe, at 
a radius of ^ 20^^ and slightly out of alignment, forming an 
axis with the core with P.A.'-^ 68°. A further, fainter exten- 
sion of the excess soft X-ray flux just reaches the position 
of the SW radio lobe at a distance of ^ AO" (^ 330 kpc) 
from the centre. The full extent of the extended emission 
along P. A.- 52° is - 70" (- 570 kpc). Unfortunately, the 
large XMM PSF and the substantial contamination of the 
soft emission by the wings of the central point source pro- 
hibit more precise characterisation of the morphology of the 
surface brightness of the extended emission. 

3.2 Spectral analysis 

We extracted X-ray spectra for the central point source and 
for the extended emission from each EPIC camera sepa- 
rately. For the core emission, we chose a circular aperture of 
10^^ radius shown in Figure lb. As the extended emission 
prohibited extraction of an annular background, suitable 
nearby regions were selected. For the extended emission, we 
extracted source and background counts from the pie-shaped 
sectors described above and shown in Figure Ic. Response 



files were created for each camera using the standard SAS 
especget tool, which we ran in extended source mode for 
the lobe emission. The spectra were grouped such that each 
bin contained 10 counts and spectra from the three cameras 
were fit simultaneously to models within the XSPEC pack- 
age. Fits were minimized using the x Gehrels statistic as 
appropriate for low count data. 90% confidence levels were 
estimated from projection of the statistic surface. A Galactic 
absorption column of Nh = 1.28 x 10^^ cm~^ was included 
in all the fits. 

A single power-law fit indicates that the unresolved ra- 
dio galaxy core is a hard X-ray source with a photon in- 
dex Tcore — l.O^Q'g. Although we cannot distinguish be- 
tween absorption and an intrinsically hard spectrum with 
the present data, such a hard spectrum usually indicates 
the presence of photoelectric absorption in an AGN. Fixing 
the photon index to a typical Tcore — 1-9 (Miyaji et al., 
2006), the best-fitting model (reduced x^=l-2) yields a sub- 
stantial column density of Nh — 1.4lj g x 10^^ cm~^ and a 
luminosity for the underlying power- lawof 1.8t°;^xl0^^ erg 
s~^, implying an aperture corrected luminosity of Lcore — 
3.6t i;^ X 10"^^ erg s ^ . The relatively poor fit statistic for the 
absorbed powerlaw results from the inability of this model 
to fit an excess of soft emission in the lowest-energy bins. 
We therefore add a second powerlaw without local obscura- 
tion to represent the soft, extended component within the 
core aperture. A much improved fit, shown in Figure |21 is 
obtained for a model with Text = 2.6. The complexity of the 
model strains the statistical limits of the dataset, however, 
and the four free parameters cannot be constrained simulta- 
neously. Freezing Text at the preferred value of 2.6 yields an 
absorption of Nh — 5.2^2 ^ ^ 10^^ cm~^, a core luminosity 
of Lcore — 3.6 X 10^^ erg s~^, and a luminosity in the soft 
power- law component of Lext — 7.3 x 10^^ erg s~^. How- 
ever, the fit can only rigourously provide a joint upper limit 
to these luminosities of 4.6 x 10^^ erg s~^. 

The statistical quality of the spectra of the extended 
source is extremely poor due to the faintness of the source 
emission. As can be seen in Figure El there are only a few 
bins in the MOS spectra and many of the bins in the pn 
spectrum are statistically equivalent to zero. The positive pn 
bins at 4-5 keV most likely result from imperfect subtraction 
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of the background. We fit a single power law model to the 
data. The best-fitting solution (reduced — 0.83) yields a 
soft photon index Fe = 3.21q g and a luminosity of 7.6 x 10^^ 
erg s~^. Exclusion of the possibly contaminated bins at 4-5 
keV does not significantly alter the fit. For completeness, 
we note that an equally good fit (reduced — 0.82) was 
obtained for a thermal Mekal model with a temperature of 
1.7 keV, a luminosity of 4.0 x 10^^ erg s~^ , and an abundance 
of 0.3 solar. The 3a limit on the luminosity of any thermal 
emission on the scale of the detected lobes outside the radio 
axis is 2.1 x lO'^^ erg s~"^. 

We can also estimate the luminosity of the extended 
component directly from the counts extracted from the pn 
camera in ^3.11 Following the procedure described in XMM 
Survey Science Centre memo SSC-LUX_TN-0059, Issue 3.0 
and adopting a F = 3 powerlaw spectrum with Galactic ab- 
sorption, we find a 0.5-2 keV pn Energy Conversion Factor 
of 6.804 X 10^^ counts cm^ erg ^. Dividing the net counts 
in the soft band by this value and the exposure time yields 
a soft-band flux of 2.6 x 10~^^ erg s~^ cm~^, implying an 
unabsorbed full-band luminosity of 7.5 it 2.0 x 10^^ erg s~^ 
which is in excellent agreement with the spectroscopically 
derived value. 

4 Discussion 

4.1 The nature of the extended X-ray emission 

The analyses above indicate soft, extended. X-ray emission 
coincident with the radio lobes of 4C 23.56 over a scale of 
~ 0.5 Mpc with a luminosity of ~ 7 x lO'^^ erg s~^. While 
there are several possible emission mechanisms for extended 
X-ray flux in HzRGs (see, e.g., Overzier et al., 2005 and 
references therein) , few can account for the morphology and 
magnitude of the flux seen in 4C 23.56. 

Extrapolating the synchrotron spectrum from the ra- 
dio waveband, we find the observed X-ray luminosity is two 
orders of magnitude higher than would be expected from 
synchrotron emission. Synchrotron self-Compton emission, 
often invoked to explain X-ray hotspots in local radio galax- 
ies (e.g., Harris et al. 2000), cannot produce X-rays over a 
large area because the synchrotron photon density can only 
exceed that of the CMB in compact regions. IC scattering 
of optical/IR photons from the central AGN (Brunetti et 
al., 2001) can not account for emission so far from the core. 
Shock heating of an ambient medium by the expanding ra- 
dio source — a proposed mechanism for the extended X-ray 
emission in PKS 1138-262 (Carilh et al., 2002) — would re- 
quire an enormous reservoir of hot, high filling- factor gas 
which is not detected in the off- axis regions. Finally, ther- 
mal emission from hot intracluster gas cannot be ruled out, 
but the surface brightness of the observed emission is much 
higher than expected and a contrived containment model 
would be needed to explain an extremely elongated mor- 
phology aligned with the radio lobes. 

We believe that in 4C 23.56, as in most previously re- 
ported HzRGs with X-ray-bright radio lobes, IC scattering 
of CMB photons is the most likely emission mechanism for 
the majority, if not the entirety, of the extended X-ray flux. 
Fourteen z > 1 HzRGs currently in the literature (see refer- 
ences in §1) have X-ray lobes on scales ranging from a few 
tens to a few hundreds of kpc with typical luminosities of 
a few times 10^^ erg s~^. The X-ray lobes in 4C 23.56 are 



the second largest discovered to datej, and the first >0.5 
Mpc-scale lobes discovered at the z > 2 epoch suspected to 
be the predominant era of massive galaxy formation. 

4.2 Energetics of the X-ray-detected lobes 

We now briefly consider the energetics of the X-ray detected 
lobes, assuming the entirety of the extended X-ray emission 
arises in IC scattering of CMB photons. Particles with a 
large Lorentz factor, 7e, will upscatter CMB photons of fre- 
quency fcMB to higher frequencies l^scat ciS 

-^7e- (1) 

l^CMB O 

Scattering to X-ray energies therefore requires electrons with 
a Lorentz factor 7e ^ 10^. If we assume that all electrons 
responsible for the X-ray lobes have this Lorentz factor, we 
can use the observed X-ray luminosity and the known CMB 
photon energy density to estimate the energy stored in the 
IC-emitting electron population. The energy flux in scat- 
tered X-rays is 



4 2 

Lx — NeUrad-'Je^^TC, 

while the energy of the relativistic electrons is 

Ee — NeJerUeC^, 



(2) 



(3) 



where Ne is the number density of relativistic electrons, ar 
is the Thompson cross-section, c is the speed of light, and rrie 
is the rest mass of the electron. Urad is the energy density of 
the CMB, which increases from the local value with redshift 
as {l-\-z)^. The energy in the scattering electron population 
is then given by 



3 LxrrieC 3 L44 

4f7rad7eCrT ~ 4 7e(l +2)4 



10*"* erg, 



(4) 



where L44 is the 2-10 keV X-ray luminosity in units of 10 
erg s~"^. 

As noted by Erlund et al. (2006), this calculation pro- 
vides only a lower limit on the energy in the scattering pop- 
ulation. Assuming a more realistic power-law distribution 
of Lorentz factors reveals dependencies of the total electron 
energy on the minimum and maximum Lorentz factors and 
on the X-ray photon index. These parameters are poorly 
constrained; adopting typical values increases the estimated 
energy by a factor of a few. In addition, this estimate con- 
siders only the energy in relativistic electrons. A substantial 
population of relativistic protons (e.g., Celotti & Fabian, 
1993; Sikora et al., 1995) could further increase the implied 
lobe energies by a factor of up to 10^. 

Converting our measured 0.5-8 keV luminosity to 2-10 
keV with F = 3, we find a lower limit to the energy stored in 
relativistic electrons in the lobes of 4C 23.56 of Ee — 8.0 x 
10^^ erg; adopting a more typical F = 2, ^ 2.2 x 10^^ erg. 
Comparing these estimates with those of Erlund et al. (2006) 
for three z 2 HzRGs we note that while the X-ray lobes 
in 4C 23.56 are at least twice as large as the lower-redshift 
sources, the minimum energy stored in IC-emitting electrons 
appears to be roughly comparable. Considering only losses 



X The giant radio galaxy 6C 0905+3955 at ^ = 1.88 exhibits 
X-ray emission associated with hotspots and lobes which are sep- 
arated by ~ 865 kpc (Blundell et al. 2006). 
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through IC emission, the radiative hfetime of the 4C 23.56 
X-ray lobes is ~ 10^ yr. 

4.3 Potential for environmental feedback 

An enormous reservoir of energy — equivalent to ~ 10^ 
supernovae — is stored in the IC-emitting electron popu- 
lation within the radio lobes of 4C 23.56. As the cooling 
timescale of this population through IC-scattering to X-ray 
wavelengths is relatively long, a large fraction of this energy 
could potentially be injected into the surroundings during a 
given epoch of AGN activity. The vast extent of the X-ray 
luminous lobes suggests a mechanism by which significant 
feedback could take place at radii of up to 250 kpc from the 
central nucleus. It is not clear, however, if the energy avail- 
able in the IC scattering electron population is sufficient to 
be a major source of feedback influencing the evolution of 
either the host galaxy or the forming cluster environment. 

Assuming a velocity dispersion of 300 km s~^, the lobe 
energy estimates derived above are 10-20 times lower than 
the binding energy of a 10^^ Mq elliptical and roughly 500 
times lower than the 5 x 10^^ erg Mq^ estimated by Benson 
et al. (2003) to be required to inhibit the excessive growth 
of such a galaxy at high redshift. In terms of heating of the 
intracluster medium, Wu, Fabian, & Nelson (2000) estimate 
non- gravitational heating of 1-2 keV per particle is needed 
to explain the observed deviation of galaxy cluster properties 
from the theoretical scaling relations. For a cluster with a 
gas mass of 10^^ M©, this is equivalent to ~ 5 x 10^^ erg 
or roughly 300 times our minimum energy estimate for the 
IC-scattering population in the lobes of 4C 23.56. 

Unless the assumption of a single Lorentz factor, je = 
10^, is dramatically wrong or there is a significant pro- 
ton component to the jet, the energy contained in the IC- 
scattering population in the lobes appears to fall short of 
the required feedback energies for the host galaxy and clus- 
ter environment by at least two orders of magnitude, even 
assuming perfectly efficient energy transfer acting over the 
required scales. In addition, the observation of significant 
cooling in IC-scattering of CMB photons is evidence that a 
significant fraction of this energy does not get injected into 
the local environment. 

Nevertheless, IC scattering could provide a very efficient 
mechanism for producing hard, ionising photons on scales far 
from the central AGN. Given reasonable assumptions for the 
halo gas density. X-ray luminosity at the level observed in 
4C 23.56 could potentially ionize a large volume of halo gas 
over scales of several hundred kpc, as outlined in Scharf et 
al. (2003). Moreover, several additional mechanisms could 
increase the feedback potential of the X-ray lobes dramati- 
cally. Energy stored within the magnetic field, found to be of 
order Ee in the Erlund et al. (2006) sources, may increase the 
total energy of the jet. More significant may be the pdV work 
associated with the expansion of lobes into the surrounding 
medium, which has been estimated to provide energies of up 
to 10^^ erg in lower redshift radio galaxies and observed to 
significantly heat the surrounding intracluster medium (e.g., 
Birzan et al., 2004; McNamara et al., 2005). Finally, bursts 
of AGN activity may occur several times over the course of a 
Gyr, re-accelerating the radio lobes and driving subsequent 
periods of environmental feedback. 
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